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The good agreement of standard solar models with helioseismology and the combined analysis of 
the solar neutrino experiments suggest that the solution to the solar neutrino problem is located in 
particle physics rather than in astrophysics. The most promising solution are neutrino oscillations, 
which usually are analyzed within the reduced 2-flavor scheme, because the solutions found therein 
reasonably well reproduce the recent data of Super-Kamiokande about the recoil-electron energy 
spectrum, zenith-angle and seasonal variations, and the event rate data of all the neutrino detectors. 
In this work, however, a survey of the complete parameter space of 3-flavor oscillations is performed. 
Basically eight new additional solutions could be identified, where the best one with Am 2 2 = 
2.7 x 10~ 10 eV 2 , Am 2 3 = 1.0 x 10" 5 eV 2 , 6i 2 = 23°, and 9i 3 = 1.3° (denoted SVO) is slightly more 
probable than any 2-flavor solution. While the 2-flavor results of the atmospheric neutrino problem 
(Am,23 > 10 -3 eV 2 ) would exclude all the 3-flavor solutions of this work, in the 3-flavor atmospheric- 
neutrino analysis the SLMA-solution with Am 2 2 = 7.9 x 10" 6 eV 2 , Am 2 3 = 2.5 x 10~ 4 eV 2 , 612 = 
1.4° , and O13 = 20° is still allowed. The relatively weak improvement of the fit using 3-flavor instead 
of 2-flavor oscillations, which appears to be due to an inconsistency of the different kind of data, 
indicates that there are possibly still systematic errors in at least one data set, or that the statistics 
is not yet sufficient. Besides, the ability of SNO and Borexino to discriminate the various 2- and 
3-flavor solutions is investigated. Only with very good statistics in these experiments the correct 
solution to the solar neutrino problem can be identified unambiguously. 

PACS numbers: 26.65.+t, 14.60. Pq, 96.60.Jw 



I. INTRODUCTION 

From the beginning of the first measurements of the 
solar neutrino flux on Earth B until the present time, 
the origin of the solar neutrino problem could not yet be 
resolved totally. While previously inaccurate or unknown 
physics used in solar-model calculations could have been 
made responsible for the discrepancy between measured 
and predicted solar neutrino flux, this kind of solution 
can presently almost be ruled out, basically of two rea- 
sons: 

Firstly, the high precision in the measurements of the 
solar p-mode frequencies and the development of helio- 
scismological inversion techniques enable the determi- 
nation of the solar sound-speed profile with high accu- 
racy 

[§• The comparison with standard solar models con- 
taining improved input physics like opacity, equation of 
state and microscopic diffusion shows an excellent agree- 
ment with the seismic models Predictions for the 
event rates in the solar neutrino detectors deduced from 
these standard solar models still are inconsistent with 
the measurements, confirming the solar neutrino prob- 
lem (Table g). 

Secondly, the three types of currently operating experi- 
ments, the chlorine detector jjj, the gallium experiments 
GALLEX/GNO § and SAGE §, and the Cerenkov- 
light counter Super-Kamiokande have different neutrino- 
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energy thresholds. This allows to determine the contri- 
bution of different parts of the solar neutrino spectrum 
to the total flux without explicitly taking into account 
solar-model calculations. From this analysis it has been 
inferred that the experimental results can be explained 
only with huge changes in the nuclear fusion rates. The 
best fit with the data is obtained even with a negative 
flux of neutrinos created in the electron-capture process 
of 7 Be ||]. Nevertheless, strong modifications of the re- 
action cross sections would be difficult to explain experi- 
mentally and theoretically. Moreover, even if a presently 
unknown physical process can account for the demanded 
changes, the resulting solar models would hardly be con- 
sistent with helioseismology Q. 

The most promising approach to the solution of the 
solar neutrino problem is an extension to the particle- 
physics standard model — neutrino mixing. Analogous 
to the CKM-mixing in the quark sector, weak and mass 
eigenstates of the neutrinos are supposed not to be iden- 
tical but connected by a unitary transformation. 

Under this assumption an initial solar electron- 
neutrino can be converted during its propagation to the 
Earth into another flavor, a /i- or r-neutrino (just-so os- 
cillations fiof ). Furthermore, the neutrinos may coher- 
ently scatter forward in solar matter (MSW-effect [pd|) 
altering the conversion probability for a certain set of 
mixing parameters. 

The possible values for the mixing parameters, with 
which the measured event rates in all detectors can be 
reproduced simultaneously, have been derived by various 
authors Jl^,[l3| , but only the oscillations between two fla- 
vors usually is taken into account. 
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Recently Super-Kamiokande has published more de- 
tailed information about the energy distribution of the 
recoiled electrons and the zenith-angle dependence of the 
neutrino-signal |IJ,|l5|. In the analysis of the 825-days 
data it became clear, that it is not possible to explain 
satisfactorily these data and the event rates of all de- 
tectors by one set of mixing parameters (see An 
excess of event rates in the high-energy bins was incon- 
sistent with the other data. 

It was the initial motivation of the present work to ex- 
amine, whether an expansion of the neutrino analysis to 
the more general 3-flavor case could resolve this discrep- 
ancy between the different types of data. However, after 
the Super-Kamiokande group reanalyzed their data and 
included new data (1117 days in total), in particular the 
excess in the energy-bin data could be diminished and 
now all kind of data can be explained simultaneously by 
2-flavor oscillations. Furthermore, the neutrino spectrum 
in the 8 B-decay has been measured recently in the labo- 
ratory |rlj| . Although within the errors the spectrum is 
in agreement with the one predicted theoretically in |L8| , 
the number of high-energetic neutrinos is overall higher 
than previously thought. Thus, the excess in the high- 
energy bins is further reduced which yields a yet slightly 
better repro ducti on of the data by the 2-flavor solutions 
(see section IV A). Nevertheless, it is presently still not 
clear, which solution to the solar neutrino problem is the 
correct one, and hence all possible solution should be de- 
duced. 

Therefore, in this work the most general case of 3-flavor 
oscillations is investigated without making any assump- 
tions about the mass scale as, for instance, inspired by 
the atmospheric neutrino problem. The latter is taken 
in various publications |20|-p2| as a constraint to investi- 
gate 3-flavor oscillations. But here, the aim is to examine 
whether the expansion to three flavors leads to new so- 
lutions with which the fits to all kind of data can be 
improved compared to the usual 2-flavor analysis. The 
implications for the atmospheric neutrino puzzle are dis- 
cussed afterwards. 

In section || the equations for 3-flavor oscillations are 
derived from the solution to the Klein-Gordon equation 
and the size of the parameter space is deduced. After 
describing the underlying solar model and the neutrino 
analysis in section III , the results for 2-flavor and 3-flavor 



oscillations are shown (section IV). Finally the ability of 
forthcoming experiments like SNO and Borexino to dis- 
criminate the various solutions are discussed (section |v|). 



A. Vacuum oscillations 

If neutrinos have mass, a mixing matrix similar to the 
Kobayashi-Maskawa matrix in the quark sector can be 
developed 
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where \v a ) denotes the weak and \vi) the mass eigen- 
states. The unitary matrix IA can be parameterized by 
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where Sy 

COS0 



are abbreviations for sin 0, 



and 

resp. cosUy (0 < 6y < tt/2) and 5 is a CP-violating 
phase pi} ], which is neglected in the following^ The 
equation of motion for a neutrino beam in vacuum obeys 
the Klein-Gordon equation for free particles (h = c = 1) 



(d t 2 -V 2 + M 2 )h(<,r)} = 0, 
where the mass matrix M 2 is defined as 
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with rrii being the mass of the neutrino mass eigenstate 
Vi. Generally the solution is given by a superposition of 
plane waves e i ( fcr_wt ' with the dispersion relation 



= k 2 



In the case of the Sun with an almost stationary neutrino 
flux, \vi(t,r)) can be expanded in components of fixed 
frequency l^iif))^ e~ lu)t . For a spherically symmetric flux 
of rclativistic neutrinos (k ~ oj) one finally gets 



(4) 



The constant "potential" u> can be removed by shifting 
the energy scale, and by using t = r Eq. (Q) can be for- 
mally written as a more familiar Schrodinger-type equa- 
tion 



II. THEORY OF NEUTRINO OSCILLATIONS 



(5) 



In the following an overview of the basic equations for 
neutrino oscillations are provided with particular empha- 
sis to the 3-flavor case. A more thorough description can 
be found e.g. in p3] or M. 



^The effect of the CP- violating phase on the analysis of neu- 
trino oscillation data has been elaborated in teal 
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with TL = Ekm ~ ^7 ■ The general solution to this equa- 



tion is 
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The probability P(r) a ^p to detect a neutrino v a with 
energy i? = u) as a neutrino of type vq at distance r = t 
from the source is therefore given by iU" 1 =U^) 
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Using Am? = m 2 — m 2 the matrix AM 2 writes as 





AM 2 = M' 2 - m\ Id 3 = | Am\ 2 

At 
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With the mass eigenstate of an electron-neutrino being 
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Eq. (M) yields for the survival probability of v^s in vac- 
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The superscript (3) denotes the case of 3-flavor mixing. 

Pe-ie depends on four quantities, two mass-squared dif- 
ferences, Am 2 2 and Am 2 3 (Am^ = Am 2 3 — Ato 2 2 ), 
and two mixing angles O12 and 613. The third mix- 
ing angle 823 does not appear in Eq. (||) and hence P^le 
is independent of this quantity. The survival probabili- 
ties for f M and v T depend on 823, but its value cannot 
be determined by solar-neutrino experiments, as in the 
energy range of the solar neutrinos and v r interact 
equally with the detector material via NC-interactions. 
Thus, only the total number of \i- plus r-neutrinos, given 
by 1—Pe-le, influences the event rates in detectors like 
Supcr-Kamiokande, SNO or Borexino. 

For oscillations between two neutrino flavors where no 
mixing into the third flavor occurs (e.g. v e «-> v^, 813 = 
0), Eq. (0) simplifies to the well known formula 
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B. Matter effect 

During the propagation of the neutrinos through the 
Sun they coherently scatter forward on the particles of 
the solar plasma. Unlike /1- and r-neutrinos, which only 
interact via NC-reactions, electron-neutrinos can addi- 
tionally couple via W-bosons to the electrons. Thus, the 
scattering cross section of a v e is altered as against the 
one of the other two neutrino-flavors. This can lead to a 
resonant flavor transition, which may create a pure v^- or 
;/ T -beam from the originally created electron-neutrinos, 
first theoretically postulated and described in |ll| (MSW- 
effeci). 

This effect can be included in Eq. (||) by substituting 
TL with TL — TL + V c ff, where 



V cS = V2G F N e rf 



influences solely the ^-contribution of the neutrino 
beam. Gf is the Fermi coupling-constant and N e the 
electron number density. The new Hamiltonian TL is no 
longer diagonal in the mass basis. To evaluate the sur- 
vival probability and thus exp(— iTLt) it is therefore nec- 
essary to diagonalize TL by a unitary transformation V 




TL,, 




v f nv. 



Similar to the vacuum oscillations the constant phase 
Mi Idz can be removed from TL as it does not change 
the survival probability. The complicated analytical ex- 
pressions for AMy have been evaluated in [^6| . Recently, 
exp(— iTLt) has been calculated in by using Cayley- 
Hamilton's theorem without explicitly deriving V. In the 
present work, however, V and Mi are computed with a 
fast numerical algorithm using none of the analytic ex- 
pressions. 



oscillations (9 defined analogously). 



C. The parameter space 

Recently various publications came up about the ac- 
tual size of the necessary parameter space covering all 
possible solutions of the solar neutrino problem [p8| j30f . 
In this section the sometimes confusing statements about 
this topic are summarized and clarified. 

The mixing angles can be defined to lie in the first 
quadrant by appropriately adjusting the neutrino field 
phases similar as in the quark sector |3l| . This can also be 

verified from the final formulae, as, for instance, P(r)e-le 
(10) depends solely on the square of sin 9^- resp. cos 

(ij — 12, 13). Moreover, it has been shown in 128] that is 

, v _ is also sufficient to consider < 923 < 17 2, if P(ry e -Lfi 

r (3) 

resp. P(r)e-> T are to be measured. In the case of the 
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matter-enhanced oscillations the situation is not as triv- 
ial, but using an analytical formula for the evolution of 
a neutrino state in matter as derived in |27j (Eq. 44 
therein), one can also show that the evolution of an ini- 
tial electron-neutrino is determined only by the squares 
of sin &ij and cos 0y- . 



1. Two flavors 

First, the case of 2-flavor oscillations is examined. Ex- 
changing the first and second row in the definition of the 
mass eigenstates (Eq. [j]) implies that Am 2 — > — Am 2 and 
sin 6 <-* cos 9 (<^> 9 — > (tt/2 — 9))[J Since the assign- 
ment of the masses ro, to the respective mass eigenstate 
Vj must not change the results, e.g. the case Am 2 > 0, 
< 9 < 7r/4 is equivalent to Am 2 < 0, 7r/4 < 9 < tt/2 
for any possible form of the Hamiltonian. Thus, without 
loss of generality 9 can be chosen to be within [0, 7r/4]. 

For pure vacuum oscillations Eq. (noh can be applied 

(2) — 

which yields that P(r) e -L e is independent of the sign 
of Am 2 and under the transformation 9 — * (tt/2 — Q). 
Hence, in this case it is even sufficient to consider 
Am 2 > 0. 

The situation is different for oscillations in matter, 
where the resonance condition for MSW-transition is 
given by 



N r , 



Am 2 /2E 
V2G F 



cos29. 
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A resonance may occur, if N les is positive, thus Am 2 > 
(9 < 7r/4). Although no resonance occurs for negative 
values of Am 2 , the matter still influences the evolution of 
the neutrino flavor composition in the solar interior (for 
a more detailed study see e.g. 26|). 

In Fig. |l| the effect of solar matter on the energy de- 
pendence of Pi-! e is demonstrated. For E/\Am 2 \ < 
10 8 MeV/eV 2 vacuum oscillations would yield an energy- 
independent survival probability for electron-neutrinos, 
as the vacuum oscillation length is small and therefore 
only a mean value of (Eq. |Io|) 
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FIG. 1. Annually averaged probability of detecting a 
neutrino produced in the solar center as u £ on Earth for 
sin 2 26 = 0.7 (0 < 9 < tt/4) with Am 2 > (solid and 
dash-dot-dot-dotted) resp. Am 2 < (dashed and dash-dotted 
line) . The two lines in each case show the probability neglect- 
ing (solid and dashed) resp. including (dash-dot-dot-dotted 
and dash-dotted) the Earth regeneration effect [BJ . 



Am 2 < is obtained for pure vacuum oscillations. How- 
ever, since vacuum oscillations in this mass range yield 
at most a suppression of the ^ e -flux of 50% (Eq. |l2|) and 
the solution to the solar neutrino problem demands, at 
least to explain the Homestake experiment, a stronger 
suppression of about 60%, for Am 2 < no reasonable 
solution can be obtained. Hence it is sufficient to con- 
sider in the 2-flavor MSW case similar to the 2-flavor 
vacuum oscillations solely Am 2 > (9 < 7r/4). 

Recently, it has been pointed out in J3(],|32| that the 
3<r-ranges of the LMA- and LOW-solutions (see below) 
extend outside this region, but additional solutions for 
the solar neutrino problems cannot be found there. In 
their analysis, however, Am 2 was fixed to be positive 
and thus < 9 < tt/2 has been examined. The region 
tt/4 < 9 < 7r/2 of their parameter space (termed "the 
dark side") is equivalent to the region Am 2 < and 
< 9 < 7r/4 discussed above. 



2. Three flavors 



would be measurable on Earth. For sin 29 = 0.7 this 

(2) 

leads to (P e _; e ) r = 0.65. In contrast to the case with 
Am 2 > (solid line), where the resonant flavor conver- 
sion diminishes the ^-contribution of the solar neutrino 
current, for Am 2 < (dashed line) the ^-portion is 
even enhanced compared to the pure vacuum-oscillation 
case. Thus, the minimum value of the z/ e -contribution for 



"''The indices 12 or 13 of the mass-squared difference 
resp. mixing angle are omitted in the 2-flavor case. 



The considerations of the 2-flavor case are now ex- 
tended to three flavors. In Fig. || the flavor space for the 
electron-neutrino is illustrated, where v e is represented 
as a (yet unknown) point on the surface of an eight of 
an unit sphere. While in the 2-flavor case the order- 
ing of the masses enables two cases to be distinguished 
(Am 2 < resp. Am 2 > 0), in the 3-flavor scenario six 
cases can be identified. But, since exchanging any axes in 
Fig. U maps the flavor space onto itself, each mass hierar- 
chy can of course be obtained from the "canonical" one 
(m 2 < m| < m 2 ) simply by exchanging the respective 
assignment of to Uj in Eq. (|l]). 
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Unlike the 2-flavor case, where for pure vacuum oscil- 
lations the parameter space could be further decreased, 
it is not possible in the general 3-flavor scenario, as 
e.g. 



'23 



'13 



transformation Am 2 2 



-A? 



'12- 




FIG. 2. Geometrical illustration of the flavor space avail- 
able for a v e in the mass basis (1/1,1/2, ^3) T . The spherical 
segment is divided into three sectors of equal size and shape, 
where each can be obtained from another by mirroring at 
the respective diagonal. Maximal mixing is given at the in- 
tersection of all three sector borders, i.e. sin 2 20i2 = 1 and 



sin 2 2ei3 = 



The cross in sector III indicates the case 



sin 20i2 = sin 20i3 = 1, the dashed border extends sector I 
to the area with 0i2,Oi3 < 7r/4. The dotted line splits 
sector I into two symmetric parts. (See text for more details.) 



and is valid for Am 2 3 w 1CU 3 eV 2 and Am 2 2 signifi- 
cantly smaller than Am 2 3 . Under these conditions the 
parameter space could be reduced considerably. 

In the present work, the most general case of 3-flavor 
oscillations of solar neutrinos is investigated, and thus 
these restrictions are not applicable here. Instead of 
performing a survey over the whole flavor space (0 < 
©12,013 < t/2) with canonical mass hierarchy, I prefer 
to consider the case Am\ 2 , Amj 3 > which covers two 
possible mass hierarchies. Thus, in this case only half 
of the total flavor space must be overviewed. In the il- 
lustration provided in Fig. || this reduced area is given 
by sector II and the lower half of sector I defined by the 
dotted line. 

To simplify the numerical survey and the analysis a 
variable transformation to the mixing angles is applied 
in each sector: Obviously Pj— ! e as defined in Eq. ([)]) is 
not symmetric under the exchange of the indices 2 <-> 3. 
However, sector I is symmetric under V2 «-> t/3, and it 
would be useful to have two quantities ©12 (©12, ©13) and 

©13(012,013) Which fulfill 

P e ^ e (Am 2 2 ,e 12 ,Am 2 3,e 13 ) = 

Pi!l e (Am 2 3 ,e 13 ,Am 2 2 ,e 12 ). 

Defining in sector I (0 < ©12 < 7r/4, tan 613 < cos ©12) 
the quantities Si2(=sin0i 2 ) and §13 by 
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§13 
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yields such a pair. With these quantities Eq. 
as =l-s?) 



(13) 



writes 



For 3-flavor neutrino oscillations in matter an exact 
analytic resonance condition can hardly be obtained be- 
cause of the complicated formula for the mass eigenval- 
ues AMij. If the masses are well separated, the 2-flavor 
resonance condition (Eq. [n]) can be applied for both 
systems by substituting the quantities (Am 2 , cos20) 
by (Am 2 2 cos©i3,cos26i 2 ) resp. (Am 2 3 , cos20i 3 ) Q. 
Hence, in this case considering ©12, ©13 < 7r/4 (sector I) 
extended to the dashed borders in sector III of Fig. ||) 
with Am 2 2 , Amf 3 > is sufficient to obtain all possible 
solutions where two resonances may occur. However, the 
solution to the solar neutrino problem may also be a com- 
bination of a non-resonant and a resonant oscillation and 
the masses can be in principle very similar, too. Thus, 
the whole parameter space must be taken into account. 

In [^0), for instance, 3-flavor oscillations were inves- 
tigated assuming canonical mass hierarchy and using a 
fixed value for Am 2 3 = 1CP 3 eV 2 in agreement with 
the atmospheric-neutrino results. They examined the re- 
maining parameter space applying an analytical formula 
for the v e survival probability, which approximates the 
solar electron-density profile by an exponential function 



P(r)B e = 1 - 
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which is unaltered under the exchange of the indices 
2^3. 

Moreover, it can be proven that the evolution matrix 
of neutrinos in matter is not changed when exchanging 
the indices 2 < — s- 3 |9| . Hence, using §12 and §13 in sector I 
yields a v e survival probability, which is constant under 
this transformation in vacuum as well as in matter. 

Similar to sector I also sector II (tt/4 < Oi 2 < 7r/2, 
tan 0i3 < sin©i 2 ) can be brought in a quadratic shape 
by using the transformation 



Sl3 



Sl3 



\A c 12 c 13 



(14) 



But this does not yield a survival probability invariant 
under the exchange of the indices 2 ^> 3, as under this 
operation sector II would be mapped on sector III. 
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III. CALCULATIONS 



A. Standard solar model 



For the following calculations my standard solar model 
( "GARSOM4" ) as described in (35) is used. It has been 
calculated using the latest input physics, equation of 
state and opacity from the OPAL-gro up [|36| , |37 ] , and nu- 
clear reaction rates as proposed by | p8| . In addition, 
microscopic diffusion of H, 3 He, 4 He, the CNO-isotopes, 
and 4 heavier elements (among them Fe) is included using 
the diffusion constants of |59| . By treating convection as 
a fast diffusive process the chemical changes due to nu- 
clear burning and diffusion (mixing) are evaluated in a 
common scheme. 

A peculiarity of GARSOM4 is the inclusion of realistic 
2d-hydrodynamical model atmospheres flcfl until an op- 
tical depth of 1000. The improvement of the high-degree 
p-mode frequencies due to a better reproduction of the 
superadiabatic layers just below the photosphere is sim- 
ilar to the one obtained by using ld-model atmospheres 
like in ||] . The advantage of using the 2d- instead of the 
ld-model atmospheres is the extension of the former to 
greater optical depths, where the stratification is already 
adiabatic and thus the solar model gets nearly indepen- 
dent of the applied convection theory. 

In Table | typical quantities of GARSOM4 using 
the convection theory developed in ]4l| ] are summa- 
rized. The predicted event rates of GARSOM4 for 
the three presently operating types of neutrino exper- 
iments GALLEX/GNO/SAGE, Homestake, and Super- 
Kamiokande are summarized in Table |H| together with 
the values obtained by Q. Since the input physics is 
very similar in both models, the predicted rates agree 
very well within the errors as quoted by Q . In addition, 
the influence of the assumed 8 B-neutrino spectrum to the 
predicted rates is shown. Using the recently measured 
spectrum JItJ yields slightly higher rates for all three 
experiments than including the theoretically predicted 
one [fl8| , which is caused by the somewhat larger number 
of high-energetic 8 B-neutrinos in the former spectrum 
and the strongly inclining detection probability toward 
higher energies. 

Figure || shows the sound-speed profile of GARSOM4 
compared to the seismic model inferred by f45}| . The 

TABLE I. Typical quantities of GARSOM4. a C m is the 
mixing-length parameter of the convection theory developed 
in EJ. Y and Z are helium and metal mass fractions, T and 
p temperature and mass density, respectively. The depth of 
the convective envelope in units of the solar radius is abbre- 
viated by -Rbcc- The indices i, s, and c denote initial, surface 
(photospheric) and central values, respectively. 
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FIG. 3. Comparison of modern solar models with the seis- 
mic model by [tt5| : shown is the relative (seismic-solar model) 
difference of sound speed for the model of reference [[l| (dashed 
line), reference @ (dash-dotted), and GARSOM4 (solid). 
The grey-shaded area indicates a conservative error range of 
seismic models according to p7[ . 



deviations of GARSOM4 from the latter are of the same 
size as standard solar models from other groups [p|fl6||. 

B. Neutrino-oscillation analysis 

Using the neutrino flux as provided by GARSOM4 the 
evolution of the initial electron neutrinos through Sun, 
space and Earth is computed taking into account oscilla- 
tions between the flavors (see Appendix). The electron- 
density profile and the radial distribution of the neutri- 
nos is taken from the solar model, too. For the electron- 
density profile of the Earth the spherically symmetric 
PREM- model |48| is applied. For each set of mixing 
parameters the neutrino energy spectrum observed on 
Earth is evaluated and folded with the detector response 



TABLE II. Expected event rates resp. neutrino flux in 
the three types of experiments GALLEX/GNO/SAGE (Ga), 
Homestake (CI), and Super-Kamiokande as predicted by the 
two solar models GARSOM4 and BP98 |L For GARSOM4 
the expected values usin g th e measured |17| (first row) and 
the theoretically derived |l8| (second row)^3-neutrino spec- 
trum are given. The last row provides the the measured values 
with their respective errors. 





Ga [SNU] 


CI [SNU] 


Super-K [lO^m^s" 1 ] 


GARSOM4 


128.7 
128.4 


7.79 
7.58 


5.18 
5.06 


BP98 


129t|? 


7 7+1.2 
'•'-1.0 


c ic + l.OO 
■->■ lo -0.72 


Experiment 


74.2 ±4.9 a 


2.56 ± 0.23 b 


2.40 ± 0.08 c 



0.975 0.275 0.020 0.245 0.018 1.57 



152 



0.713 0.188 



a Reference 
b Reference 
c Reference 
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functions. The combinations of mixing parameters which 
reproduce the measured data are found by applying a x 2 - 
analysis. 

There are four contributions to the total value of 
X 2 originating from the four different available data 
sets, the event rates, the recoil-electron energy spec- 
trum, the zenith-angle distribution, and the annual vari- 
ation. The latter three are available only from the Super- 
Kamiokande detector, while to the first one all three 
types of neutrino experiments contribute. 

For the event-rate portion the commonly used formula 
holds 



Xr 



E 



°f,exp + a i,th 



(15) 



where i denotes GALLEX/GNO/SAGE, Homestake or 
Super-Kamiokande. er 2 xp and a 2 h are the experimental 
resp. theoretical lcr-errors (Table £|. Since the input 
physics in GARSOM4 is similar to the one used in pi, 
the theoretical errors derived therein are taken for cr^. 
]V™ P are the measured event rates, which are quoted 
together with the uncertainties in Table 0. Note, that 
the Super-Kamiokande data are usually reported as 8 B- 
neutrino flux relatively to a standard solar-model predic- 
tion. Actually this number has to be understood as an 
event-rate ratio. The total number of measured events 
are divided by the theoretically expected value (e.g. from 
GARSOM4). This ratio is then often falsely taken to be 
the suppression rate of the total 8 B-neutrino flux. How- 
ever, with the energy window of the recoiled electrons 
being between 5.5^] and 20MeV no statement about the 
total number of 8 B-neutrinos below this window is pos- 
sible. Moreover, neutrino oscillations may alter the en- 
ergy spectrum of the 8 B electron-neutrino flux and as 
the scattering cross section of the neutrinos in Super- 
Kamiokande is energy-dependent, the same number of 
event rates can be obtained with different 8 B-neutrino 
fluxes. Thus in the present anal ysis the event rate in 
Super-Kamiokande following Eq. ( |A5| ) are used and not 
the total 8 B- neutrino flux. 

The recoil-electron energy spectrum is examined by 



xl = 



18 

E 



N cxp 



a, 



■N, 



th 



(16) 



where the sum extends over all 18 energy bins (Fig. ||(a)) 
and o~i te is the quadratic sum of statistical and system- 
atic errors taken from |j44| . Since the absolute value of 
the event rate in Super-Kamiokande has already been 
used in x\ the parameter a e is introduced, by which the 



s There are already data for the energy window from 5.0 to 
5.5 MeV available, but the systematic errors are still to be 
derived. 



spectrum can be normalized adequately, independent of 
the total rates. 

The contribution of the zenith-angle dependence (6 
bins. Fig. ||(b)) and seasonal- variation data (4 bins, 
Fig. H|(c)) is defined analogously (denoted in the following 
X% resp. Xa). 



IV. RESULTS 

A. Two flavors 

In a first step the solar neutrino problem is analyzed 
taking into account oscillations only between two flavors. 
Figure |] shows the allowed oscillation parameters using 
the experimentally derived 8 B-neutrino spectrum fiT| . if 
solely the event rates of the three experiments are fitted. 
Clearly the four commonly known solution-islands (see 
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FIG. 4. Confidence regions in the Am 2 -sin 2 20-plane for 
2-flavor neutrino oscillations using only the event rates. Due 
to the different mass ranges and dependences on the mix- 
ing angle, the MSW-solutions are shown in (a), the pure 
vacuum-oscillation solutions in (b). The light shaded areas 
reflect the 95.4% C.L., while the dark shaded regions show 
the 99.7% C.L. Also drawn in (a) are the constraints from 
Homestake (solid), GALLEX/GNO/SAGE (dashed), and Su- 
per-Kamiokande (dash-dotted line). The best-fit value of each 
solution is indicated by A (Table III). 
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TABLE III. Best-fit values for oscillations between two 
neutrinos taking into account solely the event rates of the 
three detector types. In the first four rows the measured 
8 B-neutrino spectrum jl^] is included, while in the last rows 
the spectrum calculated in |H| is assumed. The total number 
of d.o.f in this analysis is 1. 





Am 2 (eV 2 ) 


sin 2 26 


2 

Xtot 


SMA 


5.7 x 10" 6 


5.9 x 10" 3 


0.16 


LMA 


7.1 x 10" 5 


0.79 


0.79 


LOW 


2.5 x 1(T 8 


1.00 


5.00 


VO 


1.0 x 10~ 10 


0.95 


0.16 


SMA 


5.4 x 1(T 6 


6.1 x 10" 3 


0.21 


LMA 


7.3 x 1(T 5 


0.78 


0.81 


LOW 


5.0 x 10~ 8 


1.00 


3.70 


VO 


8.1 x 10~ n 


0.73 


0.61 



e.g. in [ p^|p^ ]) can be identified, the small-mixing (SMA) 
and large-mixing angle (LMA), the low mass-squared dif- 
ference (LOW), and the vacuum-oscillation (VO) solu- 
tions. The x 2 -values for the best-fit parameters in these 
solutions are quoted in Table 111. Apart from LOW all 
solutions have a ratio of y 2 to the number of degrees- 
of-freedom (d.o.f.) which is less than one and therefore 
these solutions are acceptable candidates as correct solu- 
tion for the solar neutrino puzzle. 



In the last four columns in Table [II the best-fit val- 
ues using the theoretically derived 8 B-neutrino spectrum 
of reference |l^] are provided. The somewhat higher ex- 
pected event rates (Table ||) result in slightly different 
mixing parameters compared to the case using the mea- 
sured 8 B-neutrino spectrum JTt| (first four rows). While 
the x 2 -values for the SMA-, LMA-, and VO-solutions 
are marginally worse with the theoretical spectrum, the 
LOW-solution give a slightly better fit to the experi- 
ments. However the changes are in all cases relatively 
small. 

While the event-rates alone favor the SMA-solution, in- 
cluding into the analysis the recoil-electron energy spec- 
trum, the zenith-angle and annual variations recorded 
by Super-Kamiokande, yields the LMA-solution as the 



best fit (Table IV). This result could also be found with 
the previous 825-day Super-Kamiokande data (see e.g. in 
p9|). In contrast to the earlier analyses, where no set 
of parameter could be found, which reproduces all sets 
of data at the same time Q, with the 1117-days data 
such simultaneous fits can be performed. In Table IV 
the best-fit values of these solutions, the respective x 1 ^ 
contributions and the ratio of Xtot to the available d.o.f. 
are provided. 

The 1, 2, and 3cr-regions (63.7 95.4, and 99.7% C.L.) 
including all available data (Fig. ||) are much bigger than 
in the case of taking solely the event rates (Fig. |^). 
Most probably neglecting the correlations between the 
different data sets of Super-Kamiokande has caused this 
growth. More detailed data published from the Super- 



TABLE IV. Best-fit values for 2-flavor neutrino oscillations 
and the respective x 2 - va hies taking into account all available 
data sets (event rates of GALLEX/GNO/SAGE, Homestake, 
and Super-Kamiokande and zenith-angle dependence, annual 
variation, and recoil-electron energy spectrum data of the lat- 
ter). In the first four rows the measured 8 B-neutrfno spectrum 
is utilize d [ p. 7 1 , while in the latter rows the spectrum calcu- 
lated in [hq] is assumed. Xh> Xe> Xz an d \ 2 a are the individual 
contributions to 

Xtot defined by Eqs. @ and (jl|). The 
last row specifies the contribution of each x 2 -P°rtion to the 
total number of d.o.f; the total number of d.o.f in this analysis 
is 26. 





Am 2 


(cV 2 ) 


sin 2 20 


xl 


xl 


xl 


Xa 


Xtot /d.o.f. 


SMA 


5.0 x 


10"° 


3.7 x 10~ J 


6.1 


16.2 


1.8 


1.1 


0.97 


LMA 


5.3 x 


10~ 5 


0.79 


0.9 


8.6 


2.0 


1.2 


0.48 


LOW 


1.1 X 


10~ 7 


0.89 


6.0 


8.6 


2.9 


0.1 


0.68 


VO 


6.9 x 


10 -io 


0.95 


6.8 


8.0 


1.9 


1.4 


0.70 


SMA 


5.4 x 


10"° 


4.2 x 10~ J 


3.3 


13.5 


2.1 


2.0 


0.79 


LMA 


5.0 x 


10 -5 


0.88 


2.0 


10.9 


1.8 


1.7 


0.63 


LOW 


1.0 x 


io- 7 


1.00 


5.9 


10.1 


2.7 


0.7 


0.75 


VO 


8.6 x 


10 -io 


0.99 


6.4 


7.9 


2.1 


1.5 


0.68 


d.o.f. 








3 


17 


5 


3 





Kamiokande collaboration are desirable to be able to per- 
form a more accurate analysis of the solar neutrino data. 
In addition, theoretical correlations between e.g. the an- 
nual and zenith-angle data should be taken into consider- 
ation. In order to show, that the most probable solution 
regions have indeed not changed drastically compared 
with the pure-rate analysis, the 10% C.L. areas are plot- 
ted in Fig. ||, too. The same region would result as the 
la-area, if the minimum Xt Q + -value would be zero instead 



of 12.5 (=0.48 x 26, Table IV) 



Independent of whether the theoretically derived or 
measured 8 B-neutrino spectrum is used, the LMA- 
solution can reproduce each single kind of data, the event 
rates, the recoil-electron energy spectrum, the zenith- 
angle and the seasonal variation acceptably. Although, 
Xt ot /d.o.f. °f the LOW- and VO-solutions are less than 
one, these solution do not lead to an acceptable fit of the 
event-rate data. Whether the SMA-solution can already 
be ruled out by the new Super-Kamiokande data depends 
on the 8 B-neutrino spectrum included in the analysis. 
While with the theoretically derived spectrum jl8), the 
event rates can just be reproduced (Xr/ 'd.o.f. ~ 1), in- 
cluding the measured one |0 hardly yields a reasonable 
fit. 

By assuming that the Chlorine rate is due to unknown 
systematic errors 30% higher than quoted and by reduc- 
ing the reaction rate for the 7 Be-proton capture by about 
15-20%, a VO-solution could be obtained which would be 
able to reproduce the event rates and the energy spec- 
trum of Super-Kamiokande fairly well |49|] . However, 
presently there is no evidence for any hidden system- 
atic uncertainties in Homestake, which would increase 
the event rate to the 3.5cr-level. 

The parameters of the best-fit LMA- and SMA- 
solutions are weakly modified when including all the data 
in the analysis (Tables HI and IV). The recoil-electron 
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FIG. 5. Confidence regions in the Am 2 -sin 2 20-plane for 
2-flavor neutrino oscillations using all data. The differ- 
ent shaded areas reflect from light to dark, 10, 63.7, 95.4, 
and 99.7% C.L. Also drawn in (a) are the constraints from 
the recoil-electron energy spectrum (dashed), zenith-angle 
(dashed-dotted) and annual variation (solid line) of the Su- 
per-Kamiokande data. The best-fit value of this analysis is 
indicated by o (Table 0); the A show the best-fit values of 
the analysis taking only the event-rate data (Table III). 



energy spectrum mainly influences the SMA-solution, 
while the zenith-angle data causes a slight shift of the 
best- fit values of the LMA-solution (Fig. a). 

For Am 2 w 10 -7 -10 -5 eV 2 the seasonal variations are 
caused by the zenith-angle variation (on the northern 
hemisphere more night data are recorded during winter 
than during summer), and thus no additional constraints 
can be obtained from the former data set in this parame- 
ter space. For smaller Am 2 the eccentricity of the Earth 
orbit leads to a "real" annual dependent signal, which 
can be used to constrain the mixing parameters. Note, 
that in contrast to the Am 2 > 10 _7 -region the seasonal 
dependence in Am 2 < 10~ 7 is now producing tiny day- 
night variations. 

Anyway, for Am 2 < 10 -9 eV 2 deviations of less than 
2% from an annually constant neutrino flux are predicted 




Jan Feb Mar Apr May Jun 
Dec Nov Oct Sep Aug Jul 

FIG. 6. Recoil-electron energy spectrum (a), zenith-angle 
(£) variation (b), and seasonal dependence of the best-fit 
LMA- (solid line), SMA- (dashed), LOW- (dash-dotted), and 
VO- (dash-dot-dot-dotted) solutions in the 2-flavor oscilla- 
tion case including all available solar neutrino data (Ta- 
ble III, lower part) compared with the measurements of 
Super-Kamiokande (1117-days data p4|) using the measured 
8 B-neutrino spectrum [fnl (solid) resp. the theoretically de- 
rived one Jl8j (dotted bins). All bins are equally spaced in 
the respective data area apart from the last bin in the energy 
spectrum, the range of which extends from 14 to 20 MeV. The 
flux variation due to the eccentric orbit of the Earth has been 
subtracted in (c). 



for the Super-Kamiokande dataf*| which is consistent 
within the errors with the recorded value. Thus, for the 
VO-solutions only very weak constraints can be obtained 



**In the analysis of the seasonal data the neutrino signal has 
been corrected for the l/r 2 -dependence of the flux. 
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from the present seasonal variation data, while for the re- 
gion of the LOW-solution these data provide important 
information. In fact, the position of the best-fit value 
of the LOW-solution has been changed by including the 
annual- variation data in the analysis (Fig. ||). 

In the regime of the VO-solution, the recoil-electron 
energy spectrum provides very stringent constraints on 
the allowed mixing parameters, excluding great part of 
the region favored by the rates. Hence, no good so- 
lution in the VO-region could be found, which repro- 
duces the recoil-electron energy spectrum as well as the 
rates recorded in GALLEX/SAGE/GNO, Homestake, 
and Super-Kamiokande. 

Nevertheless, the LMA-solution is presently the fa- 
vored solution to the solar neutrino problem, whereas 
the earlier favored SMA-solution seem to be almost ruled 
out. But still, improved statistics in the recoil-electron 
data of Super-Kamiokande is required to identify more 
reliably the correct solution to the neutrino problem. 



B. Three flavors 

In spite of the LMA-solution being able to repro- 
duce the recoil-electron energy spectrum, the zenith- 
angle and seasonal variations, and the event rates ac- 
ceptably, the analysis is extended to all three families to 
deduce whether a better fit to the data can be achieved. 
Besides, this is the physically correct treatment, which 
contains the 2-flavor case as a limiting one. 

The electron-neutrino survival probability for 3-flavor 
neutrino oscillations is determined by four quantities 



Am 2 2 , 



sm 2 29i2, Am?,, and sin 2 20i3, where the ap- 



12, amj 3 , anu sin ZO13 
propriate pairs describe each the mixing of two flavors. 
Hence a 4-dimensional parameter space has to be ex- 
amined to deduce all possible solutions. In the Ap- 
pendix the numerical realization is described with which 
the 4-dimensional parameter survey c an b e performed 
efficiently. As worked out in section II C , all possible 
solutions for MSW-solutions are obtained by consider- 
ing sector I and II (Fig. ||) with Am 2 2 , Am 2 3 > 0. 
Usingin sector I the quantity sin 2 20 13 as defined in 
Eq. ( |l3| ) instead of sin 2 2813, allows to describe the pure 
2-neutrino ^1-^3-oscillations equivalently to the v\-vi- 
case, which has been investigated thoroughly in various 
publications [[L2ljl3]. In addition, with the survival prob- 
ability P(r)i57(Am 2 2 , sin 2 20 12 , Am 2 3 , sin 2 20 13 ) being 
symmetric in the exchange of the indices 2 and 3 unnec- 
essary computations can be avoided. Nevertheless, the 
computations are very extensive due to the 4-dimensional 
parameter space. Therefore, the grid in the 3-flavor 
oscillation-survey has to be chosen less dense than in 
the 2-flavor case, where only a 2-dimensional grid had 
to be overviewed. However, the grid must still be fine 
enough that those solutions are not missed which might 
be confined to small regions in the parameter space. Fur- 
thermore, the number of neutrino paths from the Sun to 



the detector to cover the whole year of data recording 
was reduced compared to the pure 2-flavor neutrino os- 
cillations. This leads to slightly different x 2 -vahies for 
effectively pure 2-flavor solutions, which are also found 
in the full analysis. 

Subspaces of the entire possible parameter space have 

been investigated e.g. in [|o) or J50|) . In both publications 

(3) 

an analytical expression for the survival probability P e _ » e 
derived in j24|] has been used, which is valid for large mass 
separations resp. small mixing angles. In this description 

(3) (2) 

Pe-ie is determined by two 2-flavor probabilities Pe-ie for 
each mass splitting Ato 2 2 and Am 2 3 . By approximating 

the electron-density profile in the solar interior with an 

(2) 

exponential function, Pe-le can also be evaluated analyt- 
ically @. 

The results obtained in j2(J or |)(| could always been 
reproduced in the respective mass ranges. However, in 
those investigations of the 3-flavor scenario, only the 
event-rate data were available. Using the new type of 
data (recoil-electron energy spectrum, zenith-angle and 
seasonal variations) almost all solutions found in pO 50 
are disfavored. Furthermore, new solutions are identified 
since the respective parts of the parameter space were 
not covered in the analyses therein. 

In the x 2 -analysis, applied to constrain the mixing pa- 
rameters, always the whole available experimental data 
set was used for the present study. Furthermore, solely 
the measured 8 B-neutrino spectrum Jl^| has been in- 
cluded. The Am 2 2 -Amf 3 -plane is divided according 
to the 2-flavor case in three subregions where the two 
oscillation-branches v\ <-> i>2 and v\ <-> ^3 are either both 
of matter type, or both of vacuum type, or one of matter 
and the second of vacuum nature. 



1. Oscillations in the MSW mass regime 



First sector I is investigated, where V\ <-> and V\ 
both may undergo resonant MSW-transition, i.e. 10" 



» < 



Am 2 2 ,Am 2 3 < 10~ 3 eV 2 . In Fig. [?] the projection of 
the region with % 2 <26.5 on the 3-dimcnsional subspace 
sin 2 29i2-Am 2 3 -sin 2 2O13 is shown. The allowed param- 
eters are located within the grey-shaded surroundings. 

The pure 2-flavor v\ «-> v 3 oscillations where the mixing 
into Vi is negligible are given by vertical structures. For 
instance, the LMA- and LOW-solutions are represented 
by the half-pillars at sin 2 2613 = 1 an d Am 2 3 ~ 5 x 10~ D 
and 10 _7 eV 2 . The horizontal planes provide the solu- 
tions independent of ^3 , and thus are 2-flavor oscillations 
in the 1 <-> 2 system, where the influence from Am 2 2 is 
not visible. However, the symmetry of sector I in the ex- 
change of the indices 2 and 3 implies that they must be 
equivalent to solutions independent of v-i- Therefore, the 
plane around sin 2 20i 2 ss 1 is equivalent to the pillar-like 
structures at sin 2 20 13 = 1 and can thus be identified as 
an overlap of mainly LMA- and LOW-solutions. Simi- 
larly the plane at sin 2 20 i2 sa 3 x 10~ 3 represents the 
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Sector 1 




FIG. 7. Shaded volume representing the projection of the 
region with \ 2 5: 26.5 into the 3-dimensional subspace of 
sin 2 20i2-Amf 3 -sin 2 2Oi3 for MSW-transitions in the v\ 
and v\ <^vz system (sector I). The allowed values are located 
within the light shaded surroundings. The position of the 
classical 2-flavor solutions are also shown, which may be in 
the 1 <-> 2 as well as in the 1 «-> 3 system. 

SMA-solution. The 3-flavor solutions with the smallest 
X 2 -values are sited near the intersection regions of the 
horizontal planes and vertical objects, which means that 
they are at least slightly more probable than the pure 
2-flavor solutions and involve indeed all three flavors. 
In Fig. [s] the regions with x 2 < 15.4 of the present sur- 



Sectorl 




FIG. 8. Zoom into the shaded volume of Fig. [] show- 
ing the position of the SLMA-, DSMA-, DLMA-, and 
SLOW-solutions. The shaded bodies enclose the regions 
with x 2 < 15.4. The circle in the sin 2 2©i2-Am 2 3 -plane at 
sin 2 2G>i3 = 1 shows the position of the tiny solution area of 
SLOW. 



vey are shown. Due to the symmetry in v 2 v% and the 
overlap of solution islands in this projection the 6 basic 
regions in this figure belong only to four distinct solu- 
tions. Two of them can be identified as "double" SMA- 
resp. LMA-solutions, i.e. the same kind of solution in 
V\ <-> i> 2 and V\ <-> u 3 (DSMA resp. DLMA). In addi- 
tion, two solutions are combinations of a SMA-solution 
in v 1 <h> v 2 (sin 2 29i2 < 1CT 3 ) and a LMA resp. LOW 
in v x <-> V3 (sin 2 26i 3 > 10 _1 and Am 2 3 w 1CT 4 
resp. lCT 7 eV 2 .), therefore denoted SLMA and SLOW. 
The extension of the SLOW-solution in the sin 2 20 12 - 
Am 2 3 -plane at sin 2 2O13 ~ 1 is very small and thus its 
position is marked by a small circle. The assignment of 
v\ <-> u 2 to SMA and v\ <-> z/3 to LMA resp. LOW is 
ambiguous and could also be chosen vice versa. This 
is reflected by the second appearance of the SLOW- 
and SLMA-solutions in Fig. || at sin 2 2©i 2 > 10 _1 and 
sin 2 2O13 ^ 10~ 3 . The respective solution islands merge 
in this projection. 

Scclur I 




FIG. 9. Projection of the DSMA-, DLMA-, SLMA-, 
and SLOW-solutions onto the Am 12 -sin 2 2O12- or the 
Am 2 3 -sin 2 20i3-plane. The dark and bright shaded ar- 
eas show the regions of y 2 < 15.6 and 13.9. Following 
the assignment of Table M, the \ 2 5i 13.9-region of the 
DSMA-, SLMA-, and SLOW-solutions merge in the projec- 
tion onto the Am 2 2 -sin 2 2©i2-plane, but are well separated 
in the Am 2 3 -sin 2 20i3-plane (see also Fig. ^|). The region 
with sin 2 2O13 < 10 -4 has been scanned separately only 
near the DSMA-solution (1CT 6 < Am 2 2 < 10~ 4 eV 2 and 
1(T 4 < sin 2 29i2 < KT 2 ). 

In Fig. § the projection of the DSMA-, DLMA-, 
SLMA-, and SLOW-solutions into the typical 2-flavor 
planes are shown. One of the areas labelled, for in- 
stance, SLMA has to be identified with oscillations in 
V\ <-* i>i and the other with ^1-1/3 -mixing. The oscilla- 
tion parameters are provided in Table ^ together with 
the x 2 -values. Since two more parameters arc adjusted 
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in the 3-neutrino as compared to the 2-neutrino case, the 
number of degrees-of-freedom reduces from 26 to 24. 

The mixing angle of the second small mixing-angle 
branch in the DSMA-solution at sin 2 2813 is about two 
order of magnitude smaller than the usual 2-flavor SMA- 
solution (compare Figs. |j] and ^). This second SMA- 
branch causes an additional resonance at E as lOMeV 
(see Fig. [h]), which enables a better fit to the recoil- 
electron energy spectrum and the event rates (compare 
Tables IV and |y|) . Thus a solution almost as good as the 
LMA-solution could be obtained (x 2 ot /d.o.f. = 0.5). 

The DLMA- and SLMA-solutions have almost the 
same values for x%, xl±Xz an d Xa as the LMA-solution 



(cp. Tables IV and |V|). Since thus no improvement 
in explaining the solar neutrino measurements could be 
achieved by combining the LMA-solution with an ad- 
ditional oscillation in the v\ «-> ^-system, the 2-flavor 
LMA-solution itself is a favored solution in the 3-flavor 
analysis, too 

In the SLOW-solution the advantages of the 2-flavor 



SMA- and LOW-solutions are united. While with the 
SMA-solution the rates can be reproduced very well (see 
Table p^ ), the LOW-solution provides good fits to the 
energy spectrum, the annual and the zenith angle varia- 
tions (Table |Fv|). 

The best-fit value of the SLOW-solution for sin 2 26 13 
is equal to 1, and thus at the border of sector I to 
sector III. Exchanging the assignment of the indices 2 
and 3 transfers the SLOW-solution to the border be- 
tween sector I and II. In the latter sector, where 7r/4 < 
sinOi2 < 7r/2 and tan ©13 < sinOi2 a survey has been 
performed, too. However no new solutions could be iden- 
tified; solely, the "foothills" of the SLOW-solution into 
this section have been found. Hence, in the mass regime 
with 1CT 9 < Am 2 2,ATO 2 3 < 10" 3 eV 2 the region with 
7r/4 < sin 612 < 7r/2 does not provide new solutions, 
which is similar to the findings in the 2-flavor case, where 
also no additional solutions have been detected in this 
parameter range. 




10° 10 1 
E (MeV) 

FIG. 10. The annually averaged survival probability of 8 B 
electron-neutrinos detected on Earth for the deduced 3-flavor 
solutions. 



2. Three-flavor vacuum oscillations 

In the case when both mass-squared differences are in 
the vacuum-oscillation regime (Ato 2 2 ,Ato 2 3 < 10 _9 eV 2 ) 
two minima could be found, denoted DVO and DVO'. 
While DVO is located in sector I (Fig. 0), i.e. < Q 12 < 
7r/4, DVO' has been found in sector II (Fig. |l2|), where 
7r/4 < 0i2 < 7r/2. Thus although they seem to have very 
similar mixing parameters (Table |V|), there are really 
distinct solutions with different properties (see Fig. fio| ). 
Note that, while in sector I the assignment of the indices 
can be exchanged, this is not possible in sector II. If an 



> 



I 10 10 




0.4 
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FIG. 11. Projection of 



0.6 

the DVO-solution 



L0 



the 

Am? 2 -sin 2 2Q 12 - resp. Amf 3 -sin 2 2G-i£-plane. The x 2 -levels 
of the regions agree with those of Fig. |9| and the "+" -symbols 
denote the best-fit values. The best-fit value of the 2-flavor 
VO-solution reproducing solely the event rates is marked by 
A, while the position of the 2-flavor VO-solution in the com- 
plete analysis is shown by o. 
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TABLE V. The most favored solutions within the 3-flavor oscillation scenario and the respective x 2 -values. The second 
column indicates the sector in the flavor space, where the solution has been found (see also Fig. ^|). Amy is in units of eV 2 
and s^j abbreviates here sin 2 2<dij. In sector I sin 2 2O13 is defined by Eq , (ft3]), analogously in sector II by Eq. jl^). Xfli Xe, 
X% and Xa are the individual contributions to Xtot as defined by Eqs. (|l5|) and Qlq) . The number of d.o.f. is 24. 
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flavor vacuum-oscillation solutions, but still the event 
rates are fitted barely acceptably. The DVO- and DVO'- 
solutions show a compromise between the VO-solution 
(A in Fig. 11) obtained by fitting solely the event-rate 
data and the VO-solution of the complete analysis (o in 
same figure) . The mixing angles for the 3-flavor vacuum 
solution are only about one half of the usual 2-flavor VO- 
solution, as v e is now oscillating nearly equally strong 
into two other flavors. 



3. Mixed vacuum and MSW oscillations 
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of the regions agree with those of Fig. [ll]. 



exchange of v-i and v% in DVO' is desired, the mixing an- 
gles have to be transformed appropriately to lie finally in 
sector III. 

With the DVO- and DVO '-solutions a slightly better 
fit to the recoil-electron energy spectrum can be achieved 
compared to the DSMA-, DLMA-, and SLMA-solutions 
(Table |v|). However, because of the event rates being 
reproduced worse, the Xt t~ vaiues °f DVO and DVO' are 
almost equal to the values of DSMA, DLMA, and SLMA. 

Compared to the 2-flavor VO-solution an improved 
fit to the event rates has been obtained with the 3- 



The combination of oscillation between two neutrino 
flavors in the MSW mass-regime and an additional vac- 
uum oscillation into the third flavor is investigated com- 
pleting the mass ranges which have not been covered by 
the previous sections. This case has been examined in 
|(| for 10" 7 < Am 2 13 < lCT 4 eV 2 and 1CT 12 < Am 2 12 < 
10 _9 eV 2 . However, in that analysis only the event rates 
have been included, but not the recoil-electron energy 
spectrum, the zenith-angle dependence nor the annual 
variation data recorded by Super-Kamiokande. Further- 
more, the Earth-regeneration effect has been neglected. 
Therefore, this case is reinvestigated including all the 
available data and calculating the electron-neutrino sur- 
vival probabilities fully consistently including the Earth 
effect like in the previous sections. 

In this survey Am 2 3 is taken to be in the mass range 
of the MSW-solutions (1CT 9 < Am 2 3 < 10~ 3 eV 2 ) and 
Am 2 2 in the vacuum oscillation area (10 -12 < Am 2 2 < 
10 -9 eV 2 ). With these conventions about Am 2 3 , two pos- 
sibilities are conceivable for the vi <-> f 3 system, a reso- 
nant and a non-resonant oscillation. Since the masses are 
well separated, the condition for a resonance obtained in 
the 2-flavor case can be applied (Eq. nfl) 



Am 2 13 /2E 
V2G F 



cos26i 3 > 0. 



Thus, an MSW flavor transition is only possible, if the 
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FIG. 13. Recoil-electron energy spectrum (a), zenith-angle 
(£) variation (b), and seasonal dependence (c) of the best 
fit 3-fiavor solutions including all available solar neutrino 
data compared with the measurements of Super-Kamiokande. 
The experimental data agree with Fig. ^ using the measured 
8 B-neutrino spectrum |17[. 



mixing angle O13 is less than 7r/4. In the V\ «-> ^2 
system only pure vacuum oscillations occur, and thus 
©12 < 7r/4. Since the solution may be a combination of a 
vacuum oscillation in the 1 <-> 2 and a resonant resp. non- 
resonant oscillation in the 1 <-> 3 branch, the complete 
essential parameter space for 10 -9 < Amf 3 < 10 -3 eV 2 
and 10 -12 < Am 2 2 < 10 _9 eV 2 is therefore covered by 
©12 < and 6i 3 < tt/2. 

Basically two minima could be found in the parame- 
ter space of this subsurvey, which are combinations of 
a SMA- resp. a LMA-solution in v\ «-> (813 < 7r/4) 
and an additional VO-mixing in v\ «-> v-i-, therefore de- 
noted SVO and LVO. The SVO-solution is similar to the 
solution "B" found in pOfl . Interestingly, the vicinity of 
the SVO-solution has been explored intensively in 




and LVO-solutions 



and Am 13 -sin 2 2613-plane. The 



FIG. 14. Projections of the SVO 
onto the Am 2 2 -sin 2 2612 
X 2 -levels are chosen to be 13.9 (bright) and 15.6 (dark shaded) 
like in the previous figures. The circles mark the best-fit 
values of the SVO-solution, the triangles the one of the 
LVO-solution. (612,613 < tt/4) 



whereas the motivation was different: The mixing pa- 
rameters of this solution agree with the neutrino proper- 
ties predicted by a grand unification (GU) theory, where 
the neutrino masses are caused by the seesaw mechanism 
with the mass of the heaviest right-handed neutrino be- 
ing of the order of the GU-scale. 

The fit of the LVO-solution to the data is worse than 
the pure 2-flavor LMA-solution supporting the result of 
part 1 of this section, that the LMA-solution cannot be 
improved by including a third oscillation. However, the 
mass-squared difference of the LMA-branch in the LVO- 
as well as in the SLMA-solution (Ara 2 3 « 2 x lCT 4 eV 2 ) 
is almost one order of magnitude higher than that of 
the LMA-solution itself. The implication of this fact for 
the combined analysis of solar and atmospheric neutrino 
problem will be discussed further in the final section. 

The SVO-solution is presently the favored solution 
to the solar neutrino problem, in contrast to all other 
3-flavor solutions even better than the 2-flavor LMA- 
solution. SVO combines the merit of the VO-solution 
to explain the energy spectrum and the property of the 
SMA-solution to reproduce the event rates (Table III). 
Interestingly, this solution has also been found to be the 
best explanation of the solar neutrino problem, using the 
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older 504- or 708-days data sets of Super-Kamiokande. 



V. PROSPECTS FOR FUTURE EXPERIMENTS 

Future neutrino experiments like SNO or Borexino will 
provide further information to determine the correct so- 
lution to the solar neutrino problem. In the following the 
ability of these detectors to discriminate the various 2- 
and 3-flavor solutions will be discussed. 



A. Borexino 

In Fig. ^ the expected recoil-electron energy spectrum 
in Borexino of the 2- and 3-flavor solutions is shown. The 
shades and hatched areas indicate the uncertainties in the 
mixing parameters of the four 2-flavor solutions. These 
areas were obtained by allowing the mixing parameters 
to vary within the 10% confidence region of the Am 2 - 
sin 2 26-plane (see Fig. ||). 

While in the presently operating detectors, the 7 Be- 
neutrinos contribute, if at all, only a small portion to the 
total event rate, Borexino is going to measure basically 
these neutrinos (Fig. |l5|). In addition, the highest en- 
ergy bin of the recoil-electron spectrum (T e > 0.75 MeV) 
is sensitive to the neutrinos of the CNO-cycle and of the 
pep-reaction. Borexino may also measure the high en- 
ergy part of the pp-neutrino spectrum, as these may con- 
tribute about 9% to the total rate in the lowest energy 
bin (T e < 0.3 MeV). 




0.5 0.6 
T e (MeV) 

FIG. 15. The contribution of different neutrino types to the 
total event rate in Borexino in the case of no oscillations. The 
neutrino flux was derived from GARSOM4. The linestyles 
correspond to pp- (solid), pep- (dotted), 7 Be- (dashed), CNO- 
(dash-dot-dot-dotted), and 8 B- neutrinos (dash-dotted). 

In the case of the SMA-solution Borexino should detect 
only about 25% of the event rate as predicted by GAR- 
SOM4, while in the LMA-case about 65% are expected 
(Fig. |l^). Actually, for the SMA-solution the rate is ba- 
sically due to the interaction of the (or v T ) with the 



detector, as in this solution the 7 Be-neutrinos almost to- 
tally consist of ^-flavor. This leads to an unmistakable 
signal of the SMA compared to the LMA-solution: An 
excess in the lowest energy bin of about 30% compared 
to the other bins caused by the pp-neutrinos. With the 
resonance in the SMA-solution being at about 0.3 MeV, 
pp-neutrinos do, unlike 7 Be- or CNO-neutrinos, not un- 
dergo an MSW-transition. Thus, with good statistics in 
the low energy bins SMA could be clearly distinguished 
from LMA. 

One of the main advantages of this detector is the abil- 
ity to monitor seasonal variations in the neutrino signal 
- the vacuum survival probability is modified by the 
eccentric orbit of the Earth (Eq. |l^). Since the 7 Be- 
neutrinos are in contrast to 8 B-neutrinos emitted to 90% 
in a monoenergetic line at 0.862 MeV the signal in the 
detector is not smeared out like in Super-Kamiokande. 
A detailed analysis about the sensitivity of Borexino to 
detect vacuum oscillations is given in ]52j |. 

The sensitivity of the 7 Be-neutrino flux in the VO- 
regime is also visible in the zenith-angle dependence of 
the event rate in Borexino (thin solid line in Fig. |l7j(b)). 
In particular, the rates of those zenith angles, under 
which the Sun appears only around the winter or sum- 
mer solstice (| cos C| > 0.6), are influenced by the mod- 
ulation of the 7 Be-neutrino signal with the Sun-Earth 
distance. Note that unlike the Earth regeneration effect, 
which in the MSW mass-range can lead to a variation of 
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FIG. 16. Expected recoil-electron energy spectrum in 
Borexino for the best-fit 2- and 3-flavor solutions. The shaded 
and hatched areas have been obtained by allowing the mixing 
parameters to vary within the 10%-region of the 2-flavor LMA 
resp. SMA-solution shown in Fig. ^. The thin lines represent 
the best-fit 2-flavor solutions in the respective area (solid: 
SMA resp. LMA, dash-dotted: LOW, dashed line: VO). 
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FIG. 17. Annual and zenith-angle variations of the recoil 
electrons in Borexino expected for the best-fit 3-flavor solu- 
tions. The thick lines correspond to same solutions as in 
Fig. [l6| In the seasonal variation (b) the flux changes due 
to the eccentricity of the Earth's orbit have been subtracted. 



the event rate with the zenith angle only during the night 
(see thin dash-dotted line in Fig. 0(b) representing the 
LOW-solution), for VO-solutions a variation in the signal 
over the whole day is apparent. The magnitude of the 
zenith- angle variation however is smaller, as each zenith- 
angle bin is an average over many positions of the Earth 
in its solar orbit. Therefore, the seasonal variation of 
the event rates are more sensitive to detect VO-solutions 
as the zenith-angle dependence. For the best-fit mix- 
ing parameters of the MSW-solutions SMA and LMA no 
Earth regeneration takes place for the energy of the 7 Be- 
neutrinos, and thus no zenith-angle variation of the event 
rates in Borexino is expected for these cases. In case of 
the LOW-solution, however, day-night variations should 
be measured (thin dash-dotted line in Fig. (L7|(b)), 

In summary for the 2-flavor case, a total rate of about 
30 ± 5% compared to GARSOM4 is predicted for the 
SMA-solution, and about 65 ± 15% for the LMA-, LOW- 
, and VO-solutions. The latter three can be disentangles 
by using the zenith-angle and seasonal data. Great sea- 
sonal variations are expected for the VO-solution, strong 
zenith-angle dependence during the night for LOW, and 
a constant neutrino flux during the whole year for the 
LMA-solution, provided that the influence of the eccen- 



tric Earth orbit has been removed. Note, that also for 
the SMA-solution neither zenith-angle nor seasonal vari- 
ations are expected. 

The expected energy, zenith-angle, and seasonal de- 
pendences of the recoil electrons in Borexino for the 
best-fit values of the eight 3-flavor solutions are shown 
in Figs. H and [l7[ too. The DSMA and DLMA cause 
very similar signals as their 2-flavor counterparts SMA 
resp. LMA, and are therefore difficult to disentangle from 
the latter in any kind of the Borexino data. 

But the SLOW-solution leads to a clear signature in 
Borexino compared to the 2-flavor solutions. It has a 
strongly reduced event rate of almost 50% in the high- 
energy bins of Borexino, which are not observed in any 
other solution. The reason for this is the relatively high 
resonance energy of about lMeV (cf. Fig. |l^), which 
causes a strong suppression of electron-type CNO- and 
pep-neutrinos, but a weaker reduction of 7 Be electron- 
neutrinos. By way of contrast, the resonances of the 
SLMA- and SVO-solutions are about 0.8 MeV (Fig. [(J), 
and thus, the 7 Be electron-neutrino are diminished 
stronger than in the SLOW-solution, which according to 
Fig. [l^ yields smaller rates below 0.7 MeV. 

In addition, the SLOW-solution leads to a zenith-angle 
dependent signal during the night due to the Earth re- 
generation effect (Fig. [13(b)) caused by the LOW-branch 
of this solution. Thus, a large day-night variation mea- 
sured by Borexino would strongly favor LOW or SLOW. 
By additional data from the energy spectrum these solu- 
tion should then be distinguishable from each other. 

The mean event rates expected for the SLMA- and 
SVO-solutions are within a region, where all other so- 
lutions are disfavored. Since for the SVO-solution also 
annual variations are expected from its VO-branch (see 
Fig. [l7](a)), both solutions may be disentangled from each 
other and from all other solutions. 

The properties of the pure 3-flavor vacuum-oscillation 
solutions DVO and DVO' are very sensitive to the ex- 
act mixing parameters. Even a very small variation in 
the mixing parameters of these solutions may produce 
a recoil-electron energy spectrum or seasonal changes 
which resembles that of the other or the LVO-solution. 
Hence a discrimination between DVO, DVO', and LVO 
will almost be impossible with Borexino. 

These solutions are also difficult to distinguish from 
the 2-flavor VO-solution. However, the most probable 
value for Am 2 of the VO-solution is almost one order 
of magnitude larger than the VO-branch(es) of DVO, 
DVO' or LVO (compare Tables [fv| and |y|) yielding a more 
frequently varying annual signal for the VO-solution. 
Hence, measuring a weak annually varying signal with 
only one maximum per half of the year would at least 
provide a strong hint, that the solution to the solar neu- 
trino problem is rather a 3-flavor DVO-, DVO'- or LVO- 
solution than a pure 2-flavor VO solution. 
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B. SNO 



Recently, SNO has started to measure solar neutri- 
nos in the same energy window as Super-Kamiokande, 
however with a better energy resolution and the ability 
to discriminate the CC- and NC-events. Since the total 
number of neutrinos does not change by neutrino oscilla- 
tions (only their flavor), the number of NC-events has to 
agree with the no-oscillation scenario, i.e. the total num- 
ber of solar neutrinos independent of their flavor can be 
determined by this quantity^]. Therefore, no departures 
in the recoil-deuterium energy spectrum, zenith-angle or 
seasonal dependence from the standard (no-oscillation) 
values are expected for NC-events. The merit of SNO 
is the measurement of the ratio of CC- to NC-events, 
as with this quantity the ^-contribution to the total so- 
lar neutrino flux is connected. A measured CC/NC-ratio 
smaller than the standard one unambiguously proves that 
the solar neutrinos, initially u e 's, are oscillating into an- 
other flavor during their flight from the solar interior to 
the Earth. 

If the excess in the high-energy bins of Super- 
Kamiokande, which was apparent in the previous data 
published reappears in the continuing measure- 

ments, the CC/NC-ratio in SNO will enable to discrimi- 
nate between an explanation by neutrino oscillations and 
by an enhanced /iep-neutrino flux, which also was under 
discussion |^9|. While for the latter solution the hep- 
neutrinos would equally contribute to the CC- and NC- 
events leading to a constant CC/NC-ratio, for the oscil- 
lation scenario a similar excess like in Super-Kamiokande 
would be also apparent in the CC/NC-ratio of SNO. 

To illustrate the expectations for the recoil-deuterium 
energy spectrum and zenith-angle dependence of the CC- 
events in SNO, the first and second moments of these 
distributions are calculated (for definition see JT^]). The 
first moment provides the mean value and the second a 
measure of the width of the distribution. 

With very small annual or daily changes in the Super- 
Kamiokande detector, such variations arc neither ex- 
pected for CC-events (e-v e -scattering) in SNO. Thus, for 
the zenith-angle moments in SNO only small deviations 
from the standard value are predicted for the 2-flavor 
SMA-, LMA-, LOW-, and VO-solutions (Fig. g(b)). 
The largest deviations of about 4% may be measured 
in case of the LMA-solution being correct. 

The error bars for these solutions in Fig. [l^ were ob- 
tained in the same manner as for the recoil-electron en- 
ergy spectrum in Borexino. Since the best-fit values are 
not very well separated and also the error bars strongly 
overlap in the plane of the first two zenith-angle mo- 
ments, it seems to be very difficult to disentangle the 
2-flavor solution by forthcoming zenith-angle data of 



Provided, no mixing with sterile neutrinos occurs. 
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FIG. 18. Predicted values for the moments of the re- 
coil-deuterium energy (upper panel) and zenith-angle spec- 
trum (lower panel) for CC-events in SNO in case of the 3-fla- 
vor solutions. The filled symbols represent the 2-flavor solu- 
tions LMA (A), SMA (•), LOW (■) and VO (♦) of Fig.g. The 
error bars were obtained in the same manner as the hashed 
resp. shaded areas in Fig. 



SNO. A similar behavior is expected for the moments 
of the recoil-deuterium energy spectrum [ |53[ shown in 
Fig. |l^(a). However, since the energy resolution in SNO 
is better than in Super-Kamiokande (Table (v|) more sig- 
nificant data can be expected to disentangle the true so- 
lution to the solar neutrino problem. A more detailed 
discussion about the properties of SNO to discriminate 
the 2-flavor solutions can be found in . 

In the energy moments a large area is still allowed 
for the VO-solution. In case Borexino does not record 
any seasonal variations all solutions containing a VO- 
branch (e.g. LVO), and in particular the 2-flavor VO- 
solution would hardly be correct. Deviation of < —5% in 
the second energy moment would then favor the DSMA- 
solution. 

Nevertheless, from Fig. |lj it becomes obvious that, 
only if SNO is measuring a big deviation of the energy or 
zenith-angle moment from the standard (no-oscillation) 
value, one of these solution can be favored, as e.g. > 3% 
in the first zenith angle further supports LMA as correct 
solution to the solar neutrino problem. But, since SNO 
may well detect no such deviations, we need additional 
results from Borexino to exclude some of the presently 
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possible solutions. 



VI. DISCUSSION 

The most favored solutions to the solar neutrino prob- 
lem in the scope of 2- and 3-flavor neutrino oscillations 
have been investigated in this work. For this purpose the 
evolution of the neutrino state was followed numerically 
from the solar interior, where the neutrinos are produced 
in the fusion reactions, through Sun, space and Earth to 
the neutrino detectors without falling back to analytical 
approximations (cf. Appendix). 

Eight new solutions in the 3-flavor case could be found, 
among which SVO, SLMA, and DSMA are at least as 
probable as the presently favored 2-flavor LMA-solution. 
The DVO-, DVO'-, LVO, DLMA-, and, SLOW-solutions, 
with which the measured event rates, recoil-electron en- 
ergy spectrum, zenith-angle dependence and annual vari- 
ation are reproduced slightly worse than with the other 
three 3-flavor solutions, are nevertheless still good can- 
didates as solution to the solar neutrino problem. The 
new experiments Borexino and SNO, which recently have 
started to operate, and improved statistics in Super- 
Kamiokande possibly will possibly enable to discrimi- 
nate the various 2- and 3-flavor solutions. Basically, 
with Borexino the vacuum-oscillation solutions will leave 
a clear footprint in the annual variation of the event 
rates, with which they can be distinguished from MSW- 
solutions. Additionally, in case of the SLOW-solutions 
the high energy bins of the recoil-electron energy spec- 
trum will be significantly depleted compared to the other 
energy bins and a day-night asymmetry of about 10% is 
expected, which only can be reproduced by the 2-flavor 
LOW-solution. 

In this work it was investigated whether the 2-flavor so- 
lutions of the solar neutrino problem can be improved by 
including all three flavors. But the 3-flavor case cannot 
be examined totally independent of the results obtained 
for the atmospheric neutrinos. The most promising solu- 
tion to the atmospheric neutrino problem is presently the 
oscillation between and v T with the mass-squared dif- 
ference Am| 3 being approximately 10~ 3 eV 2 |5f|. This 
demands that the sum or difference of Am 2 2 and Am^ 
of the 3-flavor solutions for the solar neutrino problem is 
about 10 _3 eV 2 . However, none of the presented 3-flavor 
solutions fulfills this condition, which implies that cither 
3-flavor solutions are excluded favoring the 2-flavor LMA- 
solution (0i3 = 0) or that the origin of the atmospheric 
neutrino problem are not i/^-^-oscillations. 

But analyzing the atmospheric neutrino deficit using 
3-flavor neutrino oscillations it has been found that with 
Am| 3 w 10~ 4 eV 2 solutions exist, if Ato 2 2 and sin 2 20i2 
are in the vicinity of the SMA-solution and 4° < O13 < 
22° (Fig. 2(d) in H|). Indeed, the SLMA-solution ful- 
fills this condition as its SMA- branch (y\ ^v^) is close to 
the 2-flavor SMA-solution and the best-fit parameters of 



the LMA-branch are Am 2 3 = 2.5 x 1CT 4 and 613 = 19° 
(cf. Table |v|). This solution is close to the best 3-flavor 
solution found recently in [MJ, where the 825-days data 
of Supcr-Kamiokande have been employed and where a 
priori Am 2 3 3> 10~ 4 eV 2 has been assumed. The solu- 
tion in pi]] represents not the optimal fit to the data in 
the vicinity of the SLMA-solution, as with the best-fit 
value for Am 2 3 being 2.5 x 10~ 4 eV 2 certainly the con- 
ditions assumed therein are not fulfilled. Therefore, a 
fully consistent 3-flavor analysis of the solar and atmo- 
spheric neutrino data should be performed to find the 
most probable mixing parameters of the neutrinos. 

An additional candidate, which appears to be consis- 
tent with the atmospheric neutrino results, is the LVO- 
solution, the LMA-branch of which [y\ <-> vj) has similar 
parameters as the SLMA-solution (Table M). However, 
since the oscillation length of the VO-branch in LVO is 
much larger than the Earth diameter, v\ <-* -oscillations 
have no effect on atmospheric neutrinos. Hence, the pure 
2-flavor analysis can be performed for the atmospheric 
neutrino oscillations favoring Arn| 3 = Amf 3 « 10 _3 eV 2 . 

A further possible solution to the atmospheric neutrino 
problem may be the oscillation into a sterile neutrino, 
which is neglected in this work. So, it cannot be excluded 
that 3-flavor oscillation solutions with a sterile and an 
active neutrino- branch exist, which solve the solar and 
atmospheric neutrino problem, too. 

In the 2-flavor solutions LMA, SMA, LOW, and VO 
it is implicitly assumed that the mixing into the third 
flavor is negligible. Thus, these solutions are not cou- 
pled to the atmospheric neutrino problem, as they can 
easily be taken as z^-z^-oscillations, where no mixing 
between 24 and v r appears. If however the results of 
the LSND-collaboration JstJ, which claimed to have de- 
tected ^e-z/^-oscillations with Am 2 2 > 10 -2 eV 2 , are con- 
firmed by KARMEN2 |§], it will be difficult to solve 
the solar and atmospheric neutrino anomaly simultane- 
ously by the oscillation of three neutrino flavors. A 4- 
neutrino scheme has been proposed (one sterile and three 
active neutrinos) to cover the three different mass scales 
involved by the LSND-result, atmospheric and solar neu- 
trino problems [ p9[ . In any case, identifying the solution 
to the solar neutrino puzzle by future data of forthcom- 
ing and presently operating solar-neutrino detectors may 
also constrain the solutions for the atmospheric neutrino 
problem or vice versa. 

The possibility that 2-flavor neutrino oscillations are 
responsible for the solar neutrino puzzle appears to be 
supported by the relatively weak improvement of the fit 
by increasing the parameter space from 2 (two flavors) 
to 4 dimensions (three flavors). Indeed, only with the 
SVO-solution a somewhat better fit to the data could be 
achieved compared to the LMA-solution. This can be 
taken as a hint, that the errors in the data of Super- 
Kamiokande are greater then assumed, perhaps still un- 
known systematic errors are apparent, or that the statis- 
tics is still insufficient. Moreover, correlations between 
the different data sets of Super-Kamiokande probably are 
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important for a more accurate and reliable analysis of the 
solar neutrino data. 

Anyway, if the explanation of the atmospheric neu- 
trino deficit is, despite the presently strong evidence, 
found not to be neutrino oscillations, 3-flavor oscillations 
are allowed in any mass range, implying that the SVO- 
solution leads to the best fit to the data. Remarkably, 
SVO has been found to be the favored solution also when 
taking the previous 504- or 708-days data sets of Super- 
Kamiokande Q . Hence, SVO is together with the 2-flavor 
LMA-solution one of the stablest explanation for the so- 
lar neutrino problem. 
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APPENDIX: FROM THE INITIAL NEUTRINO 
STATE TO THE EVENT RATES 

To obtain event rates for the various detectors a va- 
riety of averages and integrals have to be evaluated. In 
this section briefly a numerical scheme is described, how 
accurate predictions for the event rates can be calculated 
numerically. 




FIG. 19. Illustrating the definition of ro and i9o- The light 
shaded region reflects the region of neutrino production, the 
dark shaded is covered by about 50 different paths (see text) 
with direction to the Earth (arrow labeled v). 



For a certain set of parameters Am 2 2 , Am 2 3 , sin 2 29i2, 
and sin 2 20 13 the evolution of an initially pure electron- 
neutrino from the solar interior to the terrestrial detec- 
tors has to be evaluated. Therefore, 

\ Vl (r)) = e- i&r \v?) (Al) 

must be integrated (cf. Eq. - \vf) denotes the ini- 
tial electron-neutrino state given by Eq. (||). The 
evolution starts at a point inside the Sun given 
by r and i? (see Fig. [l9]). The neutrino path 
has then to be followed through Sun, space and, 
if necessary, Earth arriving at the detector at the 
day D under the zenith angle £• From the fi- 
nal state \v\) the electron- neutrino survival probability 
P e ^ e (Amj 2 /E, Am 2 13 /E : sin 2 2Q 12 , sin 2 26i 3 , r , O> D, Q 
is derived using the relation 

P_ e =|k f k e }f- (A2) 

Since none of the detectors is able to discriminate be- 
tween and v T , the contribution of these neutrino fam- 
ilies is simply given by 1 — P e ->e- 

First, the neutrino paths through the Sun are consid- 
ered. Since N e , and thus also the Hamiltonian H, is a 
function of r the integration of Eq. ( ]Al| ) is performed 
piecewise. For this purpose, the interior of a solar model 
is subdivided into a grid equally spaced in electron den- 
sity (about 500 grid points). By approximating the elec- 
tron den sity between two neighboring grid point linearly, 
Eq. ( |Al| ) can be solved analytically. Thereby the neu- 
trino state can be followed through the whole solar inte- 
rior. Furthermore, the accuracy is improved by increas- 
ing the grid density near the resonance (Eq. [ll]) . 

Most of the neutrinos emerge between < ro < 0.4i? Q 
and < i)q < ir. If the MSW-resonance lies within this 
region, it is crucial to carefully calculate the paths for 
neutrinos with different ro, as the more centrally pro- 
duced neutrinos may undergo a MSW-conversion, while 
the outermost are hardly influenced by the matter. Paths 
with i}q > ir/2 may even twice cross the MSW-resonance 
region. The contribution of a neutrino emerging at an 
angle $o and radius ro to the total flux originating in 
this shell is given by 

ft(tf ) = 0.5sin(2i? )di?o 

for < i^o < 7r/2. From the symmetry around the plane 
^o = tt/2 it follows that fi(i? ) = fl(n - i? )< 

Since all neutrinos cross the shell O.4i? < ro < Rq 
a set of paths through this area is established (generally 
about 50 equally spaced in sin^o) on which exp(— iHr) 
is calculated numerically by adding up the solutions in 
each region linearized in N c . Thus, a set of matrices is 
obtained ^©(sin^o) which describes the evolution from 
r s» 0.4i?Q to the solar surface. By this means the neu- 
trino state \vf (ro, i?o)) is evaluated by first integrating 
Eq. flAl| ) until r w O.4i? and then using that matrix 
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5o(sin^o)j which was calculated for a path closest to the 
actual path of \vi(ro, $o))- It turns out that calculating 
solely two paths ($o = and n) yields mostly already ac- 
ceptable results. In this case, of course, only one matrix 
<S©(0) is needed. 

Now, the state \ vf (ro, i?o)) is followed until the surface 
of the Earth using the equation for the vacuum oscilla- 
tions (Eq. [|). Since the Earth orbit is eccentric, different 
days D during a year are discriminated (about 20 from 
perihelion to aphelion) . In a final step the rotation of the 
Earth and thus the different paths of a neutrino through 
the Earth's interior during one day are considered. Here 
a procedure similar to the shell outside the fusion region 
of the Sun is applied: A set of about 30 matrices Se (sin £) 
is calculated integrating exp(— iHr) through the Earth. 



By applying finally Eq. (A2) the survival probabil- 
ity P e ^ e (ro, D, C) for a certain set of parameters 
Aml 2 /E, Aml 3 /E, sin 2 29i 2 , sin 2 20i 3 is obtained. 
This is the basic quantity from which the final event rates 
are calculated. 

Since the experiments provide data binned in cer- 
tain periods (e.g. day, night, annual seasons or a 
whole year) resp. in certain ranges of the zenith angle, 
P e ^ e (ro,'ffo,D,h) first has to be integrated using the 
zenith-angle weights Wa(Cj D) (see e.g. fL3|) for each pe- 
riod D at the detector latitude A. 

After that, P^(ro,i?o) is folded with the radial distri- 
bution $ r (r ) of the neutrino sources inside the Sun and 
the angle weights Q discriminating the different neutrino 
types (jrp, 7 Be, 8 B etc., denoted by the superscript /). 

The event rates in the detectors are obtained by finally 
f 

folding P{ with the energy spectrum of each neutrino 

type and the detector response function Si(E), and 
summing up the contribution of all neutrino types. This 
yields 



N, 



E 



PlZ t (E<)& E (E>)dE>, 



(A3) 



where i is either Ga or CI. The detector response function 
for the Homestake detector can be found e.g. in |lq] and 
for the Ga-detectors in [p0fl . 

In Super-Kamiokande, Borexino and SNO the energy 
spectrum of the recoil electrons (resp. deuterium atoms) 
can be measured at least within certain energy bins. 
Furthermore, SNO is able to discriminate CC- and NC- 
events, while Super-Kamiokande and Borexino measure 
the total number of events arisen by any neutrino flavor. 
According to |53| the response function for these type of 
experiments with the recoil-electron energy between To 
and T\ can be described by 



j dT f dT'R(T, T') 

J T J 



,.dar(E,T') 



dT' 



(A4) 



R(T, T') = 



2ttA T ' 



exp 



(T-T + S) 2 
2A^ 



the energy resolution of the detector is included. Therein 
5 accounts for a possible uncertainty in the absolute en- 
ergy calibration and At is the energy resolution width 
following the photon statistics which yields 



A 



10 



T 



lOMeV 



with Aio being the energy resolution at lOMeV. The 
parameters Aio and 5 for Super-Kamiokande, SNO and 
Borexino are summarized in Table VI. 



TABLE VI. Energy resolution parameters used in the 
present work. The values for Super-Kamiokande and SNO 
are taken from |53|. For Borexino Ai the energy resolution 
width at 1 MeV provided by L. Oberauer (private communi- 
cation) is quoted. 





Super-K 


SNO 


Borexino 


Aio 


1.6 MeV 


1.1 MeV 


47keV 


S 


10 keV 


100 keV 





The cross sections dcrr ^ ,T - ) for Super-Kamiokande 
and Borexino were calculated following the description 
in [6]] ]. T denotes either v e and (or equivalent u T ) 
for the interaction of these neutrino flavors with Super- 
Kamiokande and Borexino or CC resp. NC for the reac- 
tion type in the SNO detector. With Er(E, Tq, T\) given 
by Eq. (A4) the event rates in Borexino resp. Super- 
Kamiokande are obtained from 



N, 



+ P^S Ve (E',T ,T 1 )\ & B (E')dE', (A5) 

while the number of CC- and NC-events in SNO is defined 
analogous to the Gallium and Chlorine detector (Eq. |A3| ). 
With the threshold energy being greater than 5 MeV only 
the sum over 8 B- and Ziep-neutrinos has to be performed 
in Super-Kamiokande and SNO. 
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